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Application to the electron gas: Exchange only 



2: 



Tim Baldsiefen, 1 ' 2 F. G. Eich, 1 ' 2 and E. K. U. Gross 2 

1 Institut fur Theoretische Physik, Freie Universitat Berlin, Arnimallee 14, D-14195 Berlin, Germany 
2 Max- Planck- Institut fur Mikrostrukturphysik, Weinberg 2, D-06112 Halle, Germany 

(Dated: August 24, 2012) 

Using the newly introduced theory of finite-temperature reduced density matrix functional theory 
[1], we apply the first-order approximation to the homogeneous electron gas. We consider both 
collinear spin states as well as symmetry broken states describing planar spin spirals and investigate 
the magnetic phase diagram as well as the temperature-dependence of the single particle spectra. 

PACS numbers: 31.15.ec,31.15.E-,65.40.-b,71.10.Ca 



INTRODUCTION 

An accurate theoretical description of the homoge- 
neous electron gas (HEG) poses a central, though not 
yet completely solved, problem in quantum mechanics. 
Highly accurate Monte-Carlo calculations of the ground- 
state energy of the HEG with collinear spin-magnetism 
form the cornerstone of the highly successful local-spin- 
density approximations LSDA [2, 3] within density func- 
tional theory (DFT). At finite temperature, the equilib- 
rium thermodynamic variables are in principle accessible 
via path- integral quantum Monte-Carlo calculations [4]. 
However, the fcrmionic sign problem renders an accu- 
rate calculation of the equilibrium properties of the HEG 
at low temperatures almost impossible [5]. Accordingly, 
a wide variety of different approaches were developed. 
These include the introduction of approximate model 
interactions [6, 7], the utilization of the dielectric for- 
mulation (employing various approximations, including 
the hypernetted chain approximation [8, 9], the modified 
convolution approximation [10, 11] and the equation-of- 
motion approach of Singwi, Tosi, Land and Sjolandcr 
[12-14]), the mapping of quantum systems to classical 
systems at finite temperature [15, 16], and the utiliza- 
tion of finite-temperature many body perturbation the- 
ory (FT-MBPT), including non-diagrammatic local field 
corrections to the random phase approximation [17-23]. 

As an alternative to these approaches, in this work 
we will consider the theoretical framework of finite- 
temperature reduced density matrix functional theory 
(FT-RDMFT) [1]. This allows a variational treatment 
of the grand potential of an arbitrary quantum system, 
employing the 1-reduced density matrix (1RDM). Com- 
pared to the conceptually similar framework of FT-DFT, 
this finite-temperature version of RDMFT is capable of 
describing the interacting kinetic energy as well as the 
exchange energy excactly and only the correlation con- 
tributions to the interaction energy and to the entropy 
need to be approximated. FT-RDMFT also allows the 
calculation of the grand potential of a HEG subject to a 
nonlocal external potential. 



The general purpose of this work is twofold. Firstly, we 
introduce the general concept of FT-RDMFT on the basis 
of the first-order functional, neglecting correlation. This 
paves the way for the treatment of correlation in Part III 
of this work [24], which in the framework of FT-RDMFT 
consists of the inclusion of an additional functional. 

Secondly, as this first-order functional yields the finite- 
temperature Hartree-Fock (FT-HF) solution, we can in- 
vestigate the temperature dependence of the different 
grand potential contributions, the magnetic phase dia- 
gram, and the single-particle spectra for the HEG in FT- 
HF for both collinear as well as spin spiral configurations, 
therefore extending previous FT-HF results [25, 26]. 



THEORETICAL FRAMEWORK 
FT-RDMFT 

We will now shortly review the most important con- 
cepts from our article on the foundations of FT-RDMFT 

[!]■ 

In the present work, we treat grand canonical quantum 
ensembles, i.e. systems in contact with a particle and a 
heat bath. For such systems the main thermodynamic 
variable is the grand potential VL = E — fj,N — 1//3S. E 
describes the internal energy, N the particle number, S 
the entropy of the system, and j3 = l/(fcgT) where ks is 
Boltzmann's constant. The chemical potential /i and the 
temperature T govern the coupling to the particle and 
heat baths respectively. A general state of a quantum 
mechanical system is described by a statistical density 
operator (SDO) D. 

D = Y,m\*i)(Vi\ ,Wi> 0,5^ = 1, (1) 



where forms a basis of the Fock space under con- 

sideration. The grand potential and entropy expressed 
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as functionals of the SDO are 

Q[D] =tv{D(H - /IN -l//3\nD} (2) 
S[D] = -tr{D hit)}, (3) 

where H is the Hamiltonian of the system and N is the 
particle number operator. The SDO minimizing Q[D] is 
given by D eq = e -P(H-^N) / tr ^P(H-^N) j n the ba- 
sis of this variational principle, Mermin [27] showed that 
there is a one-to-one correspondence between an external 
local one-particle potential and the corresponding equi- 
librium density n eq rendering a finitc-tcmperature ver- 
sion of DFT possible. However, the kinetic energy as well 
as big parts of the interaction, including the exchange en- 
ergy, are not known as explicit functionals of the density 
and have to be approximated. It is known, on the other 
hand, that the kinetic energy and the exchange energy 
can be treated exactly by using the 1RDM of the system. 
The resulting theory, RDMFT, showed some success at 
zero temperature e.g. in describing the fundamental gap 
in molecules as well as solids, including transition metal 
oxides [28-30]. The density n(r) and the 1RDM 7(1", r') 
of a quantum system in a state described by D are given 
by 



7 ^(r,r') = tr{^(r')^(r)} 
n a (r) = 7 <7(7 (r, r), 



(4) 
(5) 



with {7/v(r)} being the common field operators. By con- 
struction, the 1RDM is hermitean and can be diagonal- 
ized. We have shown in Part I [1] that for some spe- 
cial cases, the 1RDM can be seperated into two distinct 
spin contributions. As the symmetries we are consider- 
ing in this work allow such a seperation, we will write the 
1RDM in spectral representation as 



la 



(6) 



where, due to Lowdin [31], we call the eigenvalues {n ?(T } 
and eigenstates {^(x)} occupation numbers (ON) and 
natural orbitals (NO) respectively. It was shown by Cole- 
man [32] that if the NOs form a basis and the ONs fulfill 
< Jiia < l,X]icr n ^ = N then the 1RDM is cnscmblc- 
iV-representable, i.e. it corresponds to a statistical den- 
sity operator of the form (1). 

An immediate consequence of Mcrmin's proof together 
with Eqs. (4) and (5) is the one-to-one correspondence 
between D eq and 7 e9 for the case of local external po- 
tentials. In the case of nonlocal external potentials one 
needs to use Gilbert's [33] theorem to show that the 
one-to-one correspondence between D eq and "f eq still pre- 
vails. Hence one can describe the equilibrium properties 
of a grand canonical ensemble by means of functionals of 
the 1RDM. Furthermore, in contrast to zero-temperature 
RDMFT there now exists a noninteracting Kohn-Sham 



(KS)-systcm following from the simple invcrtibility of the 
Fcrmi-Dirac-distribution [34] . 
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1 4- e P(e ia -n) 



£ia - M = -3 In 
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(7) 
(8) 



The existence of a Kohn-Sham system allows the con- 
struction of a perturbative expansion of the grand poten- 
tial functional. The first-order functional for the grand 
potential is then given by 



n[ 7 ] = n k [j] + Vexth] - iiN[i\ - l/psofr] 



(9) 



where the general forms of the different functional con- 
tributions arc given in Part I [1]. 

Although the perturbative method is explicitly derived 
on the premise of describing a grand canonical ensemble, 
we have shown [1] that under the assumption of being 
in the thermodynamic limit, the functionals for the ther- 
modynamic variables of a grand canonical ensemble co- 
incide with the ones describing a canonical one. In the 
following, we are therefore able to use our method of 
FT- RDMFT for the description of an electron gas in a 
canonical ensemble. The equilibrium state is then found 
by a minimization of the free energy functional ^[7] as 
defined as 



F[ 7 ]=0[7]+MiV[7]. 



HOMOGENEOUS ELECTRON GAS 



(10) 



A central model system for the theoretical description 
of many particle quantum systems is the HEG, an ex- 
tensive review of which can be found in Rcf. [35]. The 
density p of the HEG is defined by the Wigner-Seitz ra- 
dius r s , i.e. the radius of a sphere of constant density 
which contains one electron 



1/3 



(11) 



The characteristic energy of a system at density r s is the 
Fermi energy e p with the corresponding Fermi tempera- 
ture Tp 



£ F 
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(12) 
(13) 



By choosing a basis set of certain symmetry, one re- 
stricts the domain of minimization of the free energy 
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functional to states respecting this symmetry of the sys- 
tem. In the following we will focus on two different sym- 
metries, the first one describing collincar spins and the 
second one describing a chiral spin symmetry. 



Collinear spins 

The NOs describing a system with collinear spin sym- 
metry are plane waves. Assuming furthermore that the 
charge distribution is uniform, i.e. that we are not in the 
Wigncr-crystal phase, the 1RDM becomes 



d 3 k 



(2tt) 3 

The polarization is then defined as 



ik-(r-r') 



Art + atI ' 



(14) 



(15) 



with N a = (2tt)~ 3 / d 3 kn a (k). The functional for the 
kinetic, exchange, and entropic contributions to the free 
energy are then given as 



fife [7] 

n*[7] 
Sob] 



■-E 



d 3 k /ix fc 2 

(2^ (k) T 



(16) 



J (2tt) 3 



n CT (k)?v(k') 



4- 



;k-k') 2 



(17) 



d 3 k 



r (k) ln(n CT (k)) 



(27t) 3 

+ (l-7 V (k))ln(l-7 V (k))) (18) 



In the numerical treatment, we assume the ONs to be 
constant in small volumes Vi around the k-points k^. Eqs. 
(16) - (18) then transform into sums and we arrive at the 
final expression for the free energy functional: 

F[{n la }} = n k [{n la }] + n x [{n ilT }} - 1/ pS Q [{n ia }\. 

(19) 

The kinetic energy fifc, the exchange energy Q, x , and the 
noninteracting entropy So are given by 

n k [{n la }]=J2 n ^ (2°) 

i .a 

^x[{n ia }} = -77 E n i<T n jcr K it j (21) 
£o[{n ia -}] = (n lcr ln(n la ) 

+ (l-rw)hi(l-rw))u; i , (22) 




FIG. 1. Free energy w.r.t. polarization £ and Wigner-Seitz 
radius r s at T — 5000AT. The black line denotes the equilib- 
rium free energies for a fixed r 3 . At r s « 5.8a.u. an instanta- 
neous phase transition between paramagnetic and ferromag- 
netic phases takes place. An increase in r 3 then leads to a 
continuous transition back to the paramagnetic state. 



where 



u = - 



d 3 k k 2 



P J Vi (2^) 3 2 
dk 3 



dk 3 



4tt 



p J Vi (2tt)3 J v . (2tt) 3 (k k 
d 3 k 



'Y2 



P Jv, W 



(23) 
(24) 
(25) 



We are now able to investigate the magnetic phase tran- 
sitions of a HEG at finite temperature. 

As an example, we plot the free energy at T = 5000A' 
versus polarization and density of the HEG in Figure 1. 
The positions of the free energy minima for each r s are 
denoted by the black line. When reducing the density, 
i.e. increasing r s , we encounter an instantaneous phase 
transition between a paramagnetic (vanishing polariza- 
tion) and ferromagnetic (complete polarization) equilib- 
rium phase. 

As is well known, at zero temperature, the transition 
from a paramagnetic to a ferromagnetic state can be ex- 
plained on the basis of the r~ 2 and r~ x dependencies of 
kinetic and exchange energy respectively. At small r s , the 
kinetic energy, which favours a paramagnetic configura- 
tion, is dominant. Increasing r s then increases the effect 
of the exchange energy, which favours the ferromagnetic 
state. At some critical density r c the effect of the ex- 
change contribution finally overcomes the kinetic one and 
a magnetic quantum phase transition occurs. That this 
transition is instantaneous cannot be deduced from this 
simple argument and one would have to examine the ex- 
plicit form of both kinetic as well as exchange energy. To 
explain the finite-temperature behaviour, one can essen- 
tially repeat the previous considerations, however, with 



r s (a.u.) 

FIG. 2. Equilibrium polarization £ m i n of the HEG as function 
of r s for different temperatures T. Above T = 60007^ there 
is no fully polarized equilibrium state. 

the entropy contribution included. The entropy generally 
favours a completely disordered state which in our situa- 
tion means a paramagnetic one. Also, the entropy has no 
explicit density dependence. Therefore, for increasing r s 
the entropy, compared to kinetic and exchange contribu- 
tions, becomes more and more dominant. This explains 
why for increasing r s the paramagnetic state will again 
become favourable at finite temperature. These simple 
arguments explain the existence of the phase transitions 
in Figure 1. Our results show that the first-order func- 
tional of FT-RDMFT yields an instantaneous transition 
between paramagnetic and ferromagnetic configurations 
which is not the case in real systems where the quantum 
phase transition is of second order. This result has to 
be attributed to the first-order approximation as Monte- 
Carlo results for the HEG show a continuous change of 
the order parameter, i.e. the polarization, at zero tem- 
perature [3, 6]. 

Since the entropic term in the free energy has the tem- 
perature as a prefactor, one expects the ferromagnetic 
phase to vanish faster with increasing temperature. Fig- 
ure 2 shows the validity of this argument. Interestingly, 
the nature of the phase transitions does not change, even 
at high temperatures. 

In order to study the various contributions to the free 
energy seperately we have plotted kinetic energy, ex- 
change contribution and entropy in Figure 3. We see that 
the entropy and exchange contributions always show a 
monotonically decreasing behaviour w.r.t. an increase in 
the polarization £. The kinetic energy however, which is 
known at zero temperature to be monotonically increas- 
ing with £, actually becomes decreasing for high values 
of T/Tp. 

As we will see in the following, this somewhat counter- 
intuitive effect is due to the fact that the exchange con- 
tribution hinders the temperature-induced smoothening 
of the momentum distributions and stronger so for the 
ferromagnetic configuration. To elucidate this argument 
we choose a density which will yield a ferromagnetic solu- 
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FIG. 3. Clockwise from bottom left: entropy, kinetic en- 
ergy, exchange energy, and free energy at r s = 15a. u. (Tf = 
2590K) versus the polarization £ of the HEG for differ- 
ent temperatures T. The reduced temperatures are t = 

0. 0, 1.0, 1.6, 2.2. For the free energy, we plot AF(£) = F{£) - 
f(£ = 0). 

tion at zero temperature (r s — 15a. u.). Thermodynamic 
variables for both paramagnetic and ferromagnetic con- 
figurations as functions of the temperature are shown in 
Figure 4. The curves denoted by "FD" correspond to the 
FT-MBPT expressions, i.e. the FT-RDMFT functional 
applied to Fermi-Dirac momentum distributions with the 
appropriate temperature. The FT-RDMFT functional is 
then minimized to give the curves denoted by "HF" . The 
differences of the energies for ferromagnetic and param- 
agnetic configurations are then included as the "Ai_o" 
curves. 

Ax_o = fi(£ = 1) - fi(£ = 0) (26) 

We see that in the case of a noninteracting system, i.e. 
for the "FD" curves, the kinetic energies of paramagnetic 
and ferromagnetic configurations approach each other 
but do not cross. The fact that they converge for high 
temperatures can qualitatively be understood by the con- 
cept of different "effective" Fermi temperatures T F of the 
configurations. For the paramagnetic configuration at 
zero temperature the ONs occupy two Fermi spheres of 
radii k F ^ = k' F ^ — kp, one for each spin channel. Because 
in the ferromagnetic situation the ONs are restriced to 
only one spin channel there will also be only one Fermi 
sphere with increased radius k F ^ = 2 a kp (k F ^ = 0). Be- 
cause the kinetic energies are proportional to ~Y^ a {kFa) 2 
this explains the favourisation of the paramagnetic con- 
figuration at zero temperature. An increase in tempera- 
ture will now lead to a smoothening of the Fermi sphere, 

1. e. the momentum distributions, and therefore to an 
overall increase of the kinetic energy. The quickness of 
the smoothening is determined mainly by the character- 
istic energy of the system, i.e. the Fermi energy, or, 
correspondingly, the Fermi temperature. Following from 
the arguments above, the paramagnetic configuration ex- 
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FIG. 4. Clockwise from bottom left: entropy, kinetic energy, 
exchange energy, and free energy for the HF and the noninter- 
acting functionals at r s = 15a. u. (Tp = 2590A') and £ = (0, 1) 
versus the temperature T of the HEG. The black dotted and 
the grey shaded lines denote the differences between the ferro- 
and paramagnetic configurations. The free energy differences 
are scaled by 10. The behaviour of the HF-curves can be 
explained by different "effective" Fermi temperatures. 
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FIG. 5. Phase diagram of the HEG for collinear spin configu- 
ration for the first-order FT-RDMFT functional. For T < T c , 
when increasing r s , the HEG shows both instantaneous and 
continuous phase transitions. The red dashed line denotes 
2V. 



hibits a smaller Fermi temperature when compared to the 
ferromagnetic configuration. This implies that the corre- 
sponding momentum distributions are smoothened more 
quickly which in turn lead to a relative increase of the 
kinetic energy of the paramagnetic configuration. This 
effect, hower, as we can see from Figure 4, is not big 
enough to let the kinetic energy curves cross and they 
converge for T — > oo. 

The situation changes if one includes the exchange con- 
tribution and minimizes w.r.t. the momentum distribu- 
tion. The exchange contribution for itself is known to be 
minimal for a ferromagnetic configuration with a momen- 
tum distribution describing a sharp Fermi sphere. Be- 
cause it does not couple different spin channels this also 
implies that for general polarizations, the favourable con- 
figurations are the ones describing sharp Fermi spheres 
of appropriate radii. Q x therefore counteracts the effect 
of temperature which can be interpreted as an increase 
of the effective Fermi temperature Tp. As argued, this 
increase of Tp is stronger for the ferromagnetic configura- 
tion, which, in addition to our considerations of nonintcr- 
acting systems above, then leads to a cross-over of the ki- 
netic energy curves in the FT-RDMFT treatment. These 
arguments can be applied to both entropic as well as ex- 
change contributions as well, but there the ferromagnetic 
configuration exhibits a lower value for zero temperature, 
preventing a crossing of energy curves. The behaviour of 
the free energy is more complicated because the entropy 
enters negatively. It has to be pointed out, however, that 
in the case of "FD" -momentum distributions, i.e. in first 
order FT-MBPT, an increase of temperature first leads 



to an increase in the free energy before an eventually 
monotonic decrease. This is rather unphysical and can 
be appointed to the fact that the "FD" -momentum dis- 
tribution is not acquired by any sort of variational prin- 
ciple. The "HF" -momentum distribution, one the other 
hand, is explicitly determined by a minimization proce- 
dure which leads to the qualitatively correct monotonical 
decrease of the free energy with temperature. 

After these considerations, we now calculate the equi- 
librium polarization of the HEG for a wide range of den- 
sities and temperatures. The resulting magnetic phase 
diagram is shown in Figure 5. The critical Wigner-Seitz 
radius r c , marking the zero-temperature magnetic transi- 
tion between the paramagnetic and ferromagnetic phases 
has the well known value 5. 56a. it.. On increasing the 
temperature, the ferromagnetic phase gets reduced until 
after some critical point it vanishes. The correspond- 
ing temperature T c is calculated to be at about 10500 A' 
while the critical Wigner-Seitz radius becomes 7.la.u.. 
Monte-Carlo results, on the other hand, show [2] that 
r c is expected to be about 75a. u.. An investigation of 
the temperature dependence with the help of a Stoncr 
model [6] gives T c w 80A which is close to the corre- 
sponding Fermi temperature Tp(r s = 75a. it.) = 103A. 
The fact that the critical temperature in the first-order, 
i.e. exchange only, treatment of FT-RDMFT, when com- 
pared to the Fermi temperature, is qualitatively correct 
suggests that the noninteracting entropy functional de- 
scribes big parts of the interacting entropy correctly. The 
main focus in the development of more advanced func- 
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tionals for FT-RDMFT should therefore be on the repro- 
duction of the zero-temperature critical Wigncr-Seitz ra- 
dius. We therefore expect that if one uses a temperature- 
independent correlation functional which reproduces the 
zero-temperature properties of the HEG accurately, then 
the inclusion of the noninteracting entropy will yield a 
qualitatively correct phase diagram. 



Planar Spin Spirals 

Although collinear plane waves as NOs respect the 
symmetry of the system, they do not neccessarily yield 
the lowest free energy. It was shown by Overhauser in 
1962 [36] that for the electron gas at zero temperature 
in Hartree-Fock approximation, a state describing spin- 
density waves (SDW) or charge-density waves can yield 
energies below the symmetry preserving state. We will 
investigate how this changes for pure SDWs with increas- 
ing temperature of the electron gas. 

An investigation of SDWs within zero temperature 
RDMFT was done by Eich et al.[37] and we will repeat 
the most important concepts in the following. The NOs 
describing a SDW also yields a spin channel seperable 
1RDM and are given by 




cos(^) e- 1 ^/ 2 



3 iq-r/2 



_ sin e -iqT/2 
cos (^) e iq r / 2 



ikr 



(27) 
(28) 



The spin channels now do not exhibit spin-up or spin- 
down characteristics but describe spin spirals. To dis- 
tinguish the SDW index from the collinear spin index, 
we will denote it by b in the following. The free energy 
functional can now be written as 

F[{n ib ,Oi}} = 
q 2 

rubU + — - q }Xnu - n l2 ) cos(9i)Qi 

i,b i 

- ^-^2(n ib n bj )cos 2 [ * . 3 ) 



i,3,b 



y^X n ib n jb) sin 2 



2 / - 



Ki. 



l/ff (n ib \n(n ib ) + (1 - n lb ) ln(l - n ib )\uji, (29) 



,b 



where 



Qi — — / ~TZ ~k z . 



2p Jv, (2tt)= 



(30) 



U,Kij and Wi are given in Eqs. (23)-(25). The magneti- 
zation of the HEG for SDW NOs is given by 



m(r) 



^ Acos(q • r) ^ 
A sin(q • r) 
B 



(31) 



with the two amplitudes A and B 



B = - 
2 



d 3 k 



(n k i - n k2 ) sin(6k) 



(2tt) 3 

,t Tq ("kl - n k2) C0s(6k). 
(2ir) fi 



(32) 
(33) 



In this work, we will restrict ourselves to planar SDWs, 
i.e. SDWs for which the magnetization of the HEG in z- 
direction vanishes (B = 0). As argued in [7] and [37] this 
planar configuration has a lower energy than a conical 
SDW with nonvanishing z-component. Furthermore, we 
will only consider SDWs with q |j e z . These SDWs will 
be called planar spin spirals (PSS) in the following. 

The requirement of vanishing z-component can be met 
by simple constraints on the ONs and NOs. 

n (k P -k z ) b = n(k p ,k z ) b (34) 

7T 

0(fc P ,±IM) = 2 ( X T a (k P ,\k z \)) (35) 

We can now minimize the free energy functional w.r.t. 
the ONs and orbital angles. To determine if at finite tem- 
perature a PSS state exhibits a lower free energy than the 
collinear states, we have to access the numerical values 
of these configurations in our PSS calculations. Using 
A^i = J2i n u an d N% = X^ n i2, we can define a PSS- 
polarization £ PSS as 



PSS 



Ni-N a 



(36) 




By setting q = and a.k p \k z \ = 1, the PSS-NOs become 

/ 1 _l_ 1 \ „ik-r 

(37) 
(38) 

For this set of NOs, the PSS-unpolarized state £ PSS = 
corresponds to the paramagnetic state and the PSS- 
polarized one £ PSS = 1 describes the ferromagnetic so- 
lution. 

To understand the effect of an increase of q it is in- 
structive to consider a fully polarized noninteracting sys- 
tem. We show a sketch of the g-dependence of the Fermi 
surface in Figure 6. For q = 0, the ONs describe a 
Fermi sphere of radius 2"»kp around k = 0. If one in- 
creases g, one can derive from the first three terms in 
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FIG. 6. Sketch of q-dependence of Fermi surface of a fully polarized PSS state. 



Eq. (29) and the symmetry relations (34), (35) that 
the Fermi sphere will divide symmetrically along the z- 
direction. If q supercedes 2/cf, then there will be two dis- 
tinct Fermi spheres with radius fc/r, centered at k = — q/2 
and k = q/2 respectively. 

If one now includes temperature effects, then the mo- 
mentum distribution around the Fermi surface will be 
washed-out. To reproduce the paramagnetic collinear 
configuration one therefore would have to consider the 
limit q — > oo. It is therefore possible to compare the 
free energies of a ferromagnetic collinear configuration 
(q = 0,£ PSS = 1), a fully polarized PSS configuration 
of finite q and a paramagnetic collinear configuration 
(q — > oo, £ s = 1) by varying q alone. 

The dependence of the entropy and the other thermo- 
dynamic variables on q at several temperatures is shown 
in Figure 7. The entropy displays a monotonically in- 
creasing, almost linear dependence until the two Fermi 
spheres are completely seperated (at q « 2kp). This in- 
crease of the entropy in turn leads to an increase of the 
optimal q- vector (the q for which the free energy is min- 
imal). This situation is depicted for r s = 5.50 in Figure 
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q/kF q/kF 

FIG. 7. Clockwise from bottom left: entropy, kinetic energy, 
exchange energy, and free energy at r s = 5.00a.u. (2> = 
233007?) versus the PSS q-vector. The minimal free energies 
are denoted by the arrows. 



8 where we show the free energy and the amplitude of 
a PSS-state w.r.t. q for several temperatures. For the 
temperature below some critical temperature T// , the 
free energy exhibits a minimum for finite q. If the tem- 



perature is increased above T* 



pss 



instantaneously 



jumps to a value bigger than 2k p, letting the amplitude of 
the PSS vanish. The jump in q at T C PSS therefore marks 
an instantaneous phase transition of the PSS phase where 
the amplitude of the optimal PSS is approaching zero in- 
stantaneously. In his work on spin density waves (SDW) 
in an electron gas, Ovcrhauser [36] made the conjecture 
that on increasing the temperature the optimal ampli- 
tude of the SDW-state will approach zero continuously, 
giving rise to a continuous phase transition. Our results 
disprove the validity of this conjecture for PSS-states. 

We show the dependence of the amplitude of the opti- 
mal PSS w.r.t T for several r s in Figure 9. For small q the 
amplitude increases slightly when increasing the temper- 
ature which seems to be surprising at first. Temperature 
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FIG. 8. Amplitude A and free energy F for different q 
and T for r s — 5.5a.u.. The reduced temperatures are 
t = 0.244,0.249,0.254,0.260. The amplitudes change only 
slightly when increasing the temperature, but the optimal q 
increases (arrows) discontinuously, letting the optimal am- 
plitude A op t vanish. This denotes an instantaneous phase 
transition in the PSS phase. 
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FIG. 9. Amplitudes A op t of the quilibrium PSS phases for 
different r s ({4.0, 4.5, 5.0, 5.5}a.u.) and T. The corresponding 
Fermi temperatures are {36400, 28800, 23300, 19200}^. The 
solid lines are mere guides to the eye. 



is usually expected to favour states of higher disorder and 
one could therefore expect the amplitude to be reduced. 
However, in our calculations we encounter the follow- 
ing behaviour of the Fermi surface. For q/kp < 2 the 
zero-temperature Fermi surface assumes an hourglass- 
like shape as depicted in Figure 6. We find that an in- 
crease of the temperature smoothens the Fermi surface 
which in turn leads to an increase of occupation around 
the "waist" of the hourglass. The angles 9k are usually 
closer to it/ 2 for smaller values of k z which then leads to 
an increase of the PSS amplitude A via Eq. 32. We have 
to point out, however, the numerical uncertainty of these 
findings. The energy differences for slightly changed mo- 
mentum distributions are very small but the amplitude 
shows a much stronger dependence. This also explains 
the big variance of the calculated optimal amplitudes, as 
shown in Figure 9. However, we did a careful minimiza- 
tion of the functional with several k-point-mesh refine- 
ments and believe that our numerical findings are qualti- 
tatively correct. We have argued that for small temper- 
atures the amplitude of a PSS state increases. The op- 
timal q- vector, on the other hand, changes only slightly. 
Therefore, for small temperatures the optimal amplitude 
increases. It was already mentioned in [37] that at zero 
temperature, the amplitude of the PSS decreases with 
decreasing r s . Below a critical radius r^ ss — 3.5, we 
cannot resolve the amplitude of the PSS anymore. This 
reduction of the amplitude is mainly due to the fact that 
the energy difference between the paramagnetic and fer- 
romagnetic phases increases and therefore the optimal q 
approaches 2kp (see Figure 8). When considering finite 
temperatures, this leads to a decrease of the critical tem- 
perature Tf ss , as shown in Table I. 

For those r s , for which a polarized configuration is 
favourable over an unpolarized one, a formation of a PSS 
was found to increase the free energy. 



FIG. 10. Phase diagram of the HEG in first-order- /HF- FT- 
RDMFT approximation. FM: ferromagnetic collinear phase, 
PM: paramagnetic collinear phase, PP: partially polarized 
collinear phase, PSS: planar spin spiral state. The thick 
dashed red line denotes instantaneous phase transitions. The 
dotted blue line marks the Fermi temperature Tp. 



In Figure 10 we combine the results from the collinear- 
NOs and the PSS-NOs calculations to get a schematic 
diagram of the magnetic phase diagram of the HEG. 

Considering the collinear phases, the kinetic energy 
favours a paramagnetic configuration whereas the ex- 
change contribution is minimized for a ferromagnetic one. 
At zero temperature, the different r s -dependencies lead 
to a magnetic phase transition at r c . Inclusion of the en- 
tropy, which itself favours a paramagnetic configuration, 
leads to a weakening and eventually vanishing ferromag- 
netic phase with increasing temperature. This behaviour 
and the decrease of the kinetic energy w.r.t. the polar- 
ization for high temperatures could also be explained by 
the concept of "effective" Fermi temperatures. The mag- 
netic phase transition between paramagnetic and polar- 
ized phases for increasing r s was shown to be instan- 
taneous while the transition from polarized to param- 
agnetic phases for higer r s is second-order. As can be 
seen from Figure 10, the continuous phase transitions 
occur close to the Fermi temperature Tp. The PSS 
phase shows a much stronger density dependence and 
is strongest for 4 < r s < 6. In this range we also see 
a temperature-driven instantaneous phase transition be- 
tween a PSS state and the paramagnetic state. The 
phase diagram from first-order FT-RDMFT agrees qual- 
itatively with the one derived in Ref. [7] employing a 
field-theoretical approach in combination with a contact 
interaction approximation. This method also yields a 
favourable PSS phase in the paramagnetic to ferromag- 
netic phase transition at a fixed temperature, given that 
this temperature is below some critical value. This sug- 
gests that correlation functionals in FT-RDMFT should 
preserve this general property. 

This concludes our investigation of the phase diagram 
of the HEG for the first-order FT-RDMFT functional. 

We have shown before [1] that the minimization of 
this first-order functional is equivalent to a solution of 
the FT-HF equations. In the following, we will therefore 
investigate the temperature dependence of the FT-HF 
single-particle dispersion relation for collinear as well as 
for PSS configurations. 



r c (a.u.) 


4.0 


4.5 


5.0 


5.5 


Tc(K) 


2700 


3400 


4100 


5000 



TABLE I. Critical temperatures T^ ss above which no equi- 
librium PSS-phase was found. 
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FIG. 11. Hartree-Fock dispersion relation for the HEG in 
collinear spin configuration at r s — 5a. u. for different values 
of the effective temperature t = T/Tp. The black dotted line 
denotes the noninteracting relation e(k) — k 2 /2. 



Dispersion relations 

As is well known, the HF dispersion relations for 
collinear spin configuration at zero temperature can be 
calculated analytically to yield 



e a {k) 



2 



'-fd{k/k Fa ). 



(39) 



fd{x) denotes the dimension-dependent corrections. In 
three dimensions it reads 



h{x) = - 



1 1 



4.7.: 



log 



1 



(40) 



At finite temperature, because of the entropy, the min- 
imum of the free energy functional will not exhibit pinned 
states. Therefore, at the minimum of ^[{nj}] from Eq. 
(19), the following relation holds 



OF „ dfl k dSl x 
= = — - H 

drii drii dn; 



ldSp 

f3 dn. 



(41) 



With the explict form of 5o[{nj}] of Eq. (22) this yields 



(42) 



which can easily be accessed numerically. FT-RDMFT 
can therefore be used to efficiently investigate tempera- 
ture effects on the HF-cigencncrgies, e.g. the closing of 
the energy gap between spin-channels in either collinear 
or chiral spin configurations. We have chosen a density 
at which the zero-temperature groundstate is given by a 
PSS state (r s = 5a. u.) and calculated the FT-HF dis- 
persion relation for the fully polarized HEG for various 
values of t. The results are shown in Figure 11 where, 
because of the full polarization, the Fermi surface is at 



FIG. 12. Hartree-Fock dispersion relation for the PSS phase 
in the HEG at r s = 5.00a.u. (T F = 23300iY) and q = 1.6k F - 
The red medium-dahed lines (n.i.) denote the noninteracting 
dispersions e[J(fc) whereas the other lines show the dispersion 
relation for a ferromagnetic configuration for several reduced 
temperatures t. 



kp\ = 2» kp. As the exchange contribution is quadratic 
in the ONs and does not couple the different spin chan- 
nels in collinear configuration (see Eqs. (17), (21)), the 
unoccupied band shows the noninteracting dispersion re- 
lation k 2 /2. For the occupied band at zero temperature, 
we recover the exact relations to good agreement. An 
increase of the temperature then firstly lets the "valley" 
in the dispersion relation vanish and then secondly lets 
the eigenenergies approach the noninteracting relation, 
therefore closing the HF-gap. This can be understood 
from the thermal smoothening of the momentum distri- 
bution, induced by the entropy. As NOs of very different 
quantum numbers generally have only small overlap, this 
leads to a decrease of the influence of the exchange con- 
tribution compared to the kinetic energy. 

As a second choice for the NOs of the 1RDM we 
will now consider the PSS states. The single-particle 
eigenenergies can be determined in the same way as for 
the collinear case (see Eq. (42)). The groundstate at 
r s = 5a. u. is given by a fully polarized PSS configuration 
with q = 1.6fci?, which leads to the dispersion relation 
shown in Figure 12. The noninteracting PSS energies 
e°(fc) are simply given by 



1 



° b (k) = -(k 2 ±kq) 



(43) 



As a first difference to the collinear case, we see that 
the unoccupied band at zero temperature does not follow 
the noninteracting relation. This is due to the fact that 
the exchange contribution in Eq. (29) explicitly couples 
both PSS channels. As for the collinear configuration, 
the temperature first lets the "valley" in the dispersion 
relation of the occupied band disappear. With a further 
increase of the temperature, both unoccupied as well as 
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occupied bands then approach the noninteracting disper- 
sion relation. The fact that for both collinear as well as 
PSS configurations the kink in the dispersion relation dis- 
appears much faster than the overall convergence to the 
noninteraction relation can serve as an argument, why 
for intermediate temperatures an effective mass approx- 
imation might be feasible. 

CONCLUSIONS 

In Part I of this work [1] we have introduced a new the- 
oretical framework for the description of thermodynamic 
variables of quantum systems in grand canonical equi- 
librium using the 1RDM of the system. By employing 
methods from FT-MBPT this framework allows an, in 
principle exact, construction of the grand potential func- 
tional. As a first utilization of this theory we apply the 
corresponding first-order functional, whose minimization 
we have shown [1] to be equivalent to the solution of the 
finite-temperature Hartree-Fock equations, to the HEG. 
We consider both the collinear spin phase and the PSS 
phase, possessing a chiral spin symmetry. 

When compared to a first-order treatment in FT- 
MBPT, FT-RDMFT removes the unphysical increase of 
the free energy with increasing temperature. This can 
be understood from the fact that FT-MBPT yields an 
approximation to the Green's function. This only en- 
sures that approximate results are close to the real ones 
but it docs not ensure the correct behaviour with respect 
to external parameters such as the temperature. FT- 
RDMFT, on the other hand, determines the 1RDM via a 
minimization procedure, which in the case of the first or- 
der functional ensures a decrease of the free energy with 
temperature. 

The collinear phase diagram exhibits both instanta- 
neous as well as continuous phase transitions. When cal- 
culating the amplitude of the optimal PSSs we see the ex- 
pected decrease with increasing temperature. However, 
in contrast to Ovcrhauser's proposition, this decrease 
contains a nonanlytical jump in the amplitude of the 
PSS, marking an instantaneous phase transition. Sub- 
sequently, we have employed FT-RDMFT to investigate 
the temperature dependence of the Hartree-Fock disper- 
sion relations for both collinear as well as PSS phases. 
We have found that the influence of temperature reduced 
the characteristic kink in the dispersion relation rather 
quickly whereas the assumption of the limiting noninter- 
action dispersion proceeds on a much longer temperature 
scale. 

Following from our treatment of the first-order func- 
tional of FT-RDMFT, the next step will the derivation 
of approximate functionals for the correlation part of the 
grand potential. Part III [24] of this work focusses on 
this task and will present several ways to approach this 
problem. 
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